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C^ '_ tion is solved by using the Schrodinger factorization method, the theory of intertwined 
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Rodrigues formula, a generating function and A-dependent recursion relations between 

polynomials of different orders. 
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1 Introduction 

The nonlinear differential equation 

{l + Xx'^)x-{Xx)x^ + a'^x = 0, A>0, (1) 

was studied by Mathews and Lakshmanan in [1] (see also [2] ) as an example of a non-linear oscillator 
(notice a^ was written just as a in the original paper); the most remarkable property is the existence 
of solutions of the form 

X = ^sin(a;t + cp) , 

with the following additional restriction linking frequency and amplitude 



1 + AA2 • 

That is, the equation (1) represents a non-linear oscillator with periodic solutions that were qualified 
as having a "simple harmonic form". The authors also proved that (1) is obtainable from the 
Lagrangian 

which they considered as the one-dimensional analogue of the Lagrangian density 

appearing in some nonpolynomial models of quantum field theory. 

The nonlinear equation (1) is therefore an interesting example of a system with nonlinear quasi- 
harmonic oscillations. Recently, it has been proved [3] that this particular nonlinear system can be 
generalized to the two-dimensional case, and even to the n-dimcnsional case and that these higher 
dimensional systems are superintegrable with 2n — 1 quadratic constants of motion. Moreover, 
we point out that a geometric interpretation of the higher-dimensional systems was proposed in 
relation with the dynamics on spaces of constant curvature. It was also proved the existence 
of a related A-dependent isotonic oscillator and that the two-dimensional oscillator, previously 
studied in [3], admits a superintegrable modification that corresponds to the A-dependent version 
of the Smorodinski-Winternitz sytem [4, 5]. In fact this means that the deformation introduced by 
the parameter A modifies the Hamilton-Jacobi equation but preserves the existence of a multiple 
separability. 

On the other hand, Biswas et al studied in 1973 [6] the ground state as well as some excited energy 
levels of the generalized anharmonic oscillator defined by the Hamiltonian Hm = — S jdx^ -|- x^ -|- 
Ax^™", m = 2, 3, . . ., and then they proposed the use of similar techniques for the analysis of the 
Schrodinger equation involving the potential \{x'^ /{I + gx"^)). Since then, this nonpolynomial 
potential has been extensively studied by many authors [7]- [27] from different viewpoints and 
by making use of different approaches. In many cases the term x^/(l -|- gx'^) was introduced as a 
perturbation of an initial harmonic oscillator; that is, the potential to be solved was not x^ /{1+gx^) 
by itself but x^ + Ax^/(1 -|- gx^) (exact analytical solutions have only been found for certain very 
particular values of the parameters A and gr, see e.g., Refs. [12]-[24]). It is important to note that. 



in most of these cases, the derivative part of the Schrodinger equation was the standard one, that is, 
the equation arising from a Hamiltonian H with a quadratic term of the form (l/2)p^ and leading 
to a derivative term of the form — (P' /dx^ , that is 



(P 



+ x^ + A- 



X 



^r 



En^n 



dx^ (1 + gx^) 

or the corresponding two or three-dimensional versions involving the Laplace operator in IE or IE . 

The important point is that, in the Lagrangian (2), the parameter A is present not only in 
potential x"^ / {1 -\- \ x'^) but also in the kinetic term. So, this nonlinear oscillator must be considered 
as a particular case of a system with a position-dependent effective mass [28]. In fact, it is known 
the existence of physical systems for which the mass of a nonrelativistic quantum particle may vary 
with its position (they correspond, for example, to the motion in an external field in a crystal, 
electronic properties of semiconductors, liquid crystals, etc.) and in recent years special interest 
has been drawn to generalize the standard methods of solving the Schrodinger equation to this 
new class of systems [29]- [39]. For example, some techniques have been proposed to convert the 
Schrodinger equation with a position-dependent mass into a Schrodinger equation with a constant 
mass (see [39] and references therein). Nevertheless there is an important problem at the starting 
level of quantization (that is, transition from the classical system to the quantum one) since if the 
mass m becomes a spatial function, m = m{x), then the quantum version of the mass no longer 
commutes with the momentum. Therefore, different forms of presenting the kinetic term in the 
Hamiltonian H, as for example 



T= - 



1 



m{x) 



2 , 2 
p +p 



1 



mix) 



T= - 



1 



1 



■■P 



T= - 



1 



P 



m{x) 



P 



-\Jm{x) yjm{x)- 

are equivalent at the classical level but they lead to different nonequivalent Schrodinger equations. 

In fact, the quantum version of Mathews and Lakshmanan oscillator was studied in [7] by making 
use of the following Hamiltonian 



2 



\{p\ {I - gx^)] + j^ 



gx^ 



where the notation {A, 5} = AB + BA is used. More recently another different approach has 
been discussed in Ref. [40]. It was proposed a rule for the transition from the classical system to 
the quantum one that was obtained from an analysis of the geometric properties of the function 
representing the classical kinetic energy. One of the advantages of this procedure was that the 
quantum system so obtained can be studied by using a factorization method [41]- [46] and that it 
is even endowed with the property of shape invariance [47]- [58]. 

The main objective of this article is to continue with the analysis of the quantum version of 
this particular nonlinear oscillator and to present a detailed study not only of the central question 
but also of some questions related with it. The main points to be discussed in this paper can be 
summarized in the following four points: 

• Analysis of the transition from the classical A-dependent system to the quantum one. 

This is a problem of quantization of a system with position-dependent mass and in this 
particular case the quantization procedure to be analyzed is related with the existence, at 
the classical level, of a Killing vector and also of a A-dependent invariant measure. 



• Analysis and exact resolution of the A-dependent Schrodinger equation as a Sturm-Liouville 
problem. 

In fact, not just one but two different Sturm-Liouville problems are obtained: one for A 
positive and other for A negative. This also means a discussion of the properties of the 
energy levels and of wave functions depending of the sign of A. 

• Existence of a family of A-dependent orthogonal polynomials. 

The A-dependent system we study can be considered as a deformation of the harmonic os- 
cillator in the sense that for A ^ all the characteristics of the linear oscillator must be 
recovered. The A-dependent wave functions will appear as related with a family of polyno- 
mials that can be considered as A-deformations of the standard Hermite polynomials. This 
means that the resolution of the quantum problem leads in a natural way to the study of a 
new family of orthogonal polynomials. 

• Resolution of the A-dependent Schrodinger equation by using the Schrodinger factorization 
method as an approach. 

This means a discussion of this A-dependent quantum problem as a shape-invariant problem. 
Many properties of the A-dependent orthogonal polynomials, as the generating function and 
some recursion relations, are obtained as a byproduct of this factorization approach. 

In more detail, the plan of the article is as follows: In Sec. II we first discuss the quantization 
of the classical system and then we solve the A-dependent Schrodinger equation; we analyze the 
characteristics of the solutions and we prove that they depend on a family of polynomials related 
with the Hermite polynomials; in the last subsection we prove the orthogonality of these new 
polynomials. Sec. Ill is devoted to review the main characteristics of the Schrodinger factorization 
formalism and in Sec. IV we apply this method to the A-dependent nonlinear oscillator and we 
obtain, using the operators A and A'^ , the energies En and the wave functions ^„. In Sec. V 
we study the main properties of the A-dependent "deformed Hermite" polynomials; in particular 
we prove the existence of a A-dependent Rodrigues formula, of a generating function and of A- 
dependent recursion relations between polynomials of different orders. Finally, in Sec. VI we 
discuss the results and make some final comments. 



2 A-dependent Schrodinger equation 
2.1 Quantization and Schrodinger equation 

In the general case of a system with position-dependent effective mass the transition from the 
classical system to the quantum one is a difficult problem because of the ambiguities in the order 
of the factors. Nevertheless, in the particular case of the Lagrangian (2), it was proved in Ref. 
[40] that the structure of the kinetic term suggests us a very clear and direct procedure for the 
quantization. 



Let us begin by considering the one-dimensional free-particle motion characterized by the A- 
dependent kinetic term T(A) as a Lagrangian 

1 / v^ \ 
L{x,v^;X) = Tx{x,v^) = - ( -^^^^2 ) ' (3) 

and the following nonlinear equation 

{1 + Xx^)x- {Xx)x^ = 0, 

where we point out that we admit A can take both positive and negative values; of course it is 
clear that for A < 0, A = — |A|, the function (and the associated dynamics) will have a singularity 
at 1 — |A| x^ = and we shall restrict the study of the dynamics to the interior of the interval 
x^ < 1/|A| where the kinetic energy function Tx is positive definite. 

The function T\ is invariant under the action of the vector field X{X) given by 

d 



Xx{x) = Vl + Xx 



2 

dx ' 



in the sense that we have 

Xi[Tx)=0, 

where X^^ denotes the natural lift to the velocity phase space IRxlR (the tangent bundle in differ- 
ential geometric terms) of the vector field Xx , 

d / XxVr \ d 



X\ = VT+A^ TT- + 



dx ^yJl + Xx"^^ dvx' 

In differential geometric terms this property means that the vector field Xx is a Killing vector of 
the one-dimensional metric 

ds\ = { r — K ) dx . 

^ Vl + Ax2/ 

Moving to the quantum setting, we see that this property suggests the idea of working with functions 
and linear operators defined on the space obtained by considering the real line IR endowed with 
the measure d^x given by 

which is (up to a factor) the only measure invariant under Xx [40]. This means that the operator 
Px representing the linear momentum must be self-adjoint not in the standard space ^^(IR) but in 
the space L?'{^,d^x)- 

Let us start our study with the classical Hamiltonian of the A-dependent oscillator [3] 
1 1 2 

^ = ^2;^^^' + ^2^^(tTA^)' ^- = v/l + Ax2p,, g = ma\ (5) 

Figures I and II show the form of the potential Vx{x) for several values of A (A < in Figure I and 
A > in Figure II). 



As it has been pointed out in [40], the generator of the infinitesimal 'translation' symmetry 
\/l + Ax^ d/dx is skew-selfadjoint in the space L'^(M,,d^\) and therefore the transition from the 
classical system to the quantum one is given by defining the operator 



P. 



d 



so that 



ifi vT+Ax^-— , 
dx 



d 



{1 + \x^)pI ^ 
in such a way that the quantum version of the Hamiltonian (5) becomes 



d 



n' [ViTx^ ^){Vl + X^ -^) , 



H 



(1 + Ax'') -— 7 - ( )Xx-- + (-)q. , ^ 



Let us consider the following quantum Hamiltonian 
H = 



1 + Ax ] -r-^ — ( Ax- — \- (-jmaia -\ A), , „ 

2m ^ ^ dx2 ^2m^ dx ^2^ ^ m ^Vl + Ax^ 



(6) 



where we have slightly modified the value of the parameter g that now is given hy g = ma^ + XTia; 
this is done in order that the notation used in this approach coincides with the one to be presented 
in the next section III. It is also convenient to simplify this function H by introducing adimensional 
variables (y, A) defined by 



n 



ma 



A 



ma 



A, 



in such a way that the following equality is satisfied 

l + Ax^ = 1 + Ay^ 



The Hamiltonian H takes then the following form 
H = 



1 /-, * 2^ d"^ 

■^i^ + W' 



dy' 



V 



(i)Ay-f + (i)(l + A) , 

2 dy 2 V 1 + A y^ 



(ha) 



(7) 



and then the Schrodinger equation 
reduces to the following adimensional form 



e{ha) , 



^ y ^ ' 2 4^ ^ dy ^2'^ ^Vl + Ay2 



y 



2 '' ' "" ' dy"^ 
which, after a small simplification, can be finally rewritten as follows 



^ = e^ 



d^ 



d 



(l + A!,= )^* + Ai,^*-(l + A),j^^_^, 



y 



^ + (2e)^ = 0. 



It is known that in the A = case the asymptotic behaviour at the infinity suggest a factorization 
for the wave function. The idea is that a similar procedure can be applied to this A-dependent 



6 



equation. Let us first denote by ^oo the following function ^oo = (1 + Ay^)^ i/(2A) ^j^^t vanish in 
the limit y^ ^ oo (in the case A > 0) or in the limit y^ -^ — 1/A (in the case A < 0). Then we 
obtain 



d^ . d ...../ y2 



('+^''^'£?+'^='£-<'+^'(i 



dy 



+ Ay^ 



^r 



^r 



Thus, ^oo is the exact solution in the very particular case e = 1/2 and represents the asymptotic 

behaviour of the solution in the general case. Consequently, this property suggest the following 

factorization 

^(y,A) = %,A)(l+Ay2)-i/(2A), 



and then the new function h{y,A) must satisfy the differential equation 

(1 + A y2)/i" + (A _ 2)yh' + (2e - l)/i = , h = h{y, A) . 

Let us first remark that, as limA^o(l + Ay"^) = 1 and limA^o[l/(2 A)] = oo, we have 



(9) 
(10) 



Iim(l+Ay2)i/{2A) 



exp 



lim 



log(l + Ay^ 
2A 



exp 



lim 

A->0 



yV(l + Ay2 



oVV^ 



Consequently, 



limA_o^(y,A) = %,0)e-(i/2),2_ 



Furthermore, if we put A = in the equation (8) we see that h{y) = h{y, 0) satisfies the Hermite 
equation 

A^O h" -2yh' + {2e-l)h = 

In other words, for A = 0, the function /i(y, 0) is a solution of Hermite equation. Recall that when 
2e — 1 is an even integer number, i.e. 2e — 1 = 2p, such differential equation admits polynomial 
solutions, called (with an appropriately chosen overall factor) Hermite polynomials. 

Nevertheless, it is also interesting to note that the dependence of the parameter A makes of (10) 
an equation with certain similarities with the Legendre's equation. 

It is clear that the three coefficient functions of the linear differential equation 

(1 + A y2)/i" + (A - 2)yh' + (2e - l)/i = , h = h{y, A) , 

are analytic at the origin (and the director coefficient does not vanish in a neighbourhood of y = 0). 
The origin is therefore an ordinary point and we expect an analytic solution h with a power series 
expansion convergent in an interval (— i?, R) for some non-zero value of the radius of convergence 
R, 

oo 

/i(y. A) = ^ a„(A) y" = ao(A) + ai(A) y + a2(A) y^ + . . . 

n=0 

In this way the equation becomes 

oo 

^ [(n + 2)(n+l)a„+2 + Aany" + (A-2)na„ + (2e-l)a„] y" = , (11) 

n=0 



and therefore it determines to the following A-dependent recursion relation 



«n+2 = (-1) 



n + 2)(n + 1) 



n(An- 2) + (2e - 1) 



n = 0,1,2,... 



Note that this relation shows that, as in the A = case, even power coefficients are related among 
themselves and the same is true for odd power coefficients. In both cases, having in mind that 



lim, 



an+2X 



n+2 



hm^ 



SO that the radius of convergence R is therefore given by 

1 



n{Kn- 2) + (2e - 1) 
(n + 2)(n + l) 



X 



K\\x' 



R 



vm 



Hence, when we consider the limit A ^ 0, we recover the radius i? = oo of the Hermite's equation. 

The general solution h is given by the linear combination h = a^hi +aih-2 where hi{y) and h2{y) 
are the solutions determined by /ii(0) = 1, /ii(0) = and /i2(0) = 0, ^-2(0) — 1) respectively. 

The condition for the differential equation to admit a polynomial solution is the existence of an 
integer number p such that 

2e- 1 = 2p-Ap^, 



because then 



Op/O , ap+2 = , 



and the solution is a polynomial of order p. This relation establishes some possible values for the 
spectrum of the Hamiltonian H and therefore of H (this question will be discussed in the next 
subsection). 

The polynomial solutions are given by 
• Even index (even power polynomials) 

'H2p 
0'2r 



r=p 

hp Yl ^2r y 
r=0 
Co 



2r 



= (-1)' ^P' (P - 2)(p' - 4) . . . ip' - 2(r - 1)) 

'^' [2-Ap'][2-A(p' + 2)][2-A(p' + 4)]...[2-A(/ + 2(r-l))] 

where we have introduced the notation p' = 2p. More specifically, the expressions of the first 
solution hi, in the particular cases of p' = 0, 2, 4, 6, are given by: 

TCo = ko 

W2 = k2[l-2{l-A)y^] 

W4 = fe4[l-4(l-2A)y2 + (l)(i_2A)(l-3A)y4] 

He = A;6[l-6(l-3A)y2 + 4(l-3A)(l-4A)/ + (-^)(l-3A)(l-4A)(l-5A)y6] 

15 



where ko, ki, k2, ■ ■ ■, are overall multiplicative constants. 



Odd index (odd power polynomials) 

a2r+l 



r=p 

r=0 
\r+l '^ 



-1 



2r + l! 



(p'-l)(/-3)...(p'-(2r-l)) 

[2-A(p' + l)][2-A(p' + 3)]...[2-A(p' + (2r-l))] 



where we have introduced the notation p' = 2p+ 1. More specifically, the expressions of the 
second solution /i2 for p' = 1, 3, 5, are given by: 



Hi 



hy 



A:3b-(-)(l-2A)y3] 

W5 = A:5[y-(^)(l-3A)|/3 + (l)(i_3A)(l-4A)y5] 

where fci, fcs, k^, . . ., are arbitrary multiplicative constants. 

It is clear that these particular polynomials solutions Tim, m = 0,1,2, ... , play, for this A-dependent 
oscillator, a similar role to the Hermite's polynomials for the standard harmonic oscillator. 



2.2 A-dependent Sturm-Liouville problem and orthogonality 

The A-dependent differential equation 

aoh" + aih' + 02/1 = , 

ao = {l + Ay^), ai = (A-2)y, 03 = (2e - 1) , 

is not self-adjoint since Og 7^ ai but it can be reduced to self-adjoint form by making use of the 
following integrating factor 

^(y) = ^L)J(ai/ao)dy ^ ^^ ^ ^y2y{A+2)/(2 A) 

in such a way that we arrive to the following expression 



d_ 
dy 



p(y,A)^l +(2e-l)r(y,A)/i = 0, 
dy i 



where the two functions p = p{y. A) and r = r(y, A) are given by 
p{y,A) = e/("i^o)d?/^(i + ^y2^i/2-i/A^ 



oo 1 + A y^ 



l+Ay2)l/2-l/A. 



Thus, we have obtained 



(it 



Vl + Ay2 s dh 



+ 



(2e- l)h 



dyl\{l + Ay^y/^J dyi ^1 + Ay^ (1 + Ay2)i/A 



(12) 



that, together with appropriate conditions for the behaviour of the solutions at the end points, 
constitute a Sturni-Liouville problem. It is to be pointed out that the boundary conditions are 
in fact different according to the sign of A; therefore we arrive to, no just one, but two different 
Sturm-Liouville problems: 

• If A is negative the range of the variable y is limited by the restriction y^ < 1/|A|. In this case 
the problem, defined in the bounded interval [— a a, a a] with a a = l/\/|A[, is singular because 
the function p{y,A) vanishes in the two end points yi = — ca and y2 = OA- The conditions 
to be imposed in this case lead to prescribe that the solutions h{y, A) of the problem must 
be bounded functions at the two end points, yi = — (ia and 2/2 = oA; of the interval. It is 
clear that this leads to the above mentioned polynomial solutions. 

• If A is positive the variable y is defined in the whole real line IR and, therefore, the Sturm- 
Liouville problem is singular. The solutions h{y,A) must be well defined in all IR, and the 
boundary conditions prescribe that the behaviour of these functions when y ^ ±00 must be 
such that their norms, determined with respect to the weight function r{y), be finite. It is 
clear that in this case the solutions of the problem are again the A-dependent polynomials 
Hm, m = 0,1,2, ... 

Proposition 1 The eigenf unctions of the problem, (12) are orthogonal with respect to the function 

r = (l + Ay2)-{i/2+i/A). 

Proof: This statement is just a consequence of the properties of the Sturm-Liouville problems. 

Because of this the polynomial solutions Tirm m- = 0, 1, 2, . . ., of the equation (10), satisfy 

^"A n„,{y,A)nn{y,A) ^ 

— ay = u, niy^n, A<0, (13) 

and 

^"^ n^{y,K)nn{y,K) ^ 

.(m?jvvTTAF''=°' "^^"' ''^'- ^''^ 

If we define the A-dependent Herniite functions ^m by 

^m{y, A) = Hmiy, A) (l + A y^)- i/(2A) ^ ^ = 0, 1, 2, . . . 

then the above statement admits the following alternative form: The A-dependent Hermite func- 
tions ^miy,A) are orthogonal with respect to the weight function r = 1/ \Jl + Ay"^: 

^„(y,A)^„(y,A)r(y,A)dy= / ^„(y, A) *„(y. A)— =^== = , m^n, A<0, 
-OA J-aA ^/l + Ay^ 

(15) 
and 

/oo dy 

*„(y,A)^„(y,A)— =^== = 0, m^n, A>0. 
-00 vl + Ay^ 

(16) 

Note that this orthogonality with respect to the weight function r coincides with the orthogonality 
with respect to the measure d^\ discussed in the subsection 2.1. 

Finally, we also remark that 

lim A^o^m(y, A) = H^{y) e~ (^Z^) y\ m = 0, 1, 2, . . . 

10 



2.3 Wave functions and energy levels 

We started by introducing the factorization (9) for the wave functions ^ and then, when the 
boundary conditions have been taken into account, we have obtained that the first factor h{y, A) 
must be a polynomial function of a very particular class. The result is that the wave functions 
become "a A-dependent polynomial divided by the function (1 + Ay^)-*^''^^^"; more specifically, 

^„,{y, A) = nUy, A) (1 + Ay2)- ^/(^a) , ^ = 0, 1, 2, . . . 

where m is the order of the polynomial. Since the second factor is even and has no nodes, the 
corresponding eigenfunction ^m(?/;A) has the same parity as m and m nodes (when A < the 
nodes are inside the the bounded interval [— aA,aA]). Therefore, these two particular properties, 
parity and number of nodes, are preserved by the deformation introduced by A. 

The two cases A negative and A positive are rather different and must be studied separately. 

(i) Let us suppose that A is negative, A < 0. 

In this case the dynamics is restricted to the interval y^ < 1/|A| and the boundary conditions 
imply that the function ^(y, A) must vanish at the end points yi = — l/-\/|A| and ?/2 = l/-\/|A|. This 
means that the first factor h{y, A) must be a polynomial and, because of this, the eigenfunctions 
are given by a countable family of functions ^m, m = 0,1,2, .. .. The eigenfunctions \l'm and the 
eigenvalues em are given by 

^„(y,A)=H„(y,A)(l-|A|/)i/(2|A|), 

m + ^^+^m^ \A\ , m = 0,l,2,...,m,... (17) 

These eigenfunctions ^^ are well defined and vanish at the end points, yi and y2, for all the 
integers values of m without any restriction; thus, there exists a infinite but countable number of 
eigenfunctions. The energy spectrum is unbounded but, contrarily to what happens in the case of 
the linear oscillator, is not equispaced 

eo < ei < 62 < 63 < . . . < e„ < e^+i < ■■■ 

Cra+i -em = 1 + (m+ -j |A| . 

Nevertheless, it is also true that the energy eo of the ground-state ^o is eo = 1/2; that is, the 
zero-point energy is given by Eq = {l/2)h a for any A 7^ 0, as in the linear case. 

(ii) Let us suppose that A is positive, A > 0. 

This case is rather delicate and deserves to be analyzed with more detail. It is convenient to 
recall that the domain of y is now the whole real line IR and hence it is necessary to take into 
account the problem of the convergence at the infinity. In fact, it is necessary that the following 
integral be convergent 

nl{y,A) 



/ r: — : av < 00 

i-oofl+Aw2)l/A,/l + Aw2 



'-oo(l+Ay2)i/A^l + Ay2 

and, as for large values of y, the powers of the dominant terms in the numerator and the denominator 
are 2m and 1 + 2/A respectively, we arrive to a certain condition to be satisfied by m. In fact, 

11 



given a certain value of A, then the admissible functions ^m are those associated to integer values 

of m satisfying the condition 

1 

m < ruA = -7- • 

This means that for large values of A the oscillator only admits the fundamental level ^0 as eigen- 
function, while for smaller values of A the oscillator admits a certain finite number of eigenstates, in 
such a way that when A decreases this number increases and in the linear limit A ^ the number 
goes to infinity. The eigenfunctions ^^ and the eigenvalues e^ are given by 

^„(y,A)=H„(y,A)(l+Ay2)-i/(2A)^ 

m+^)-^m^A, m = 0,l,2,...,NA, (18) 

where N^ denotes the greatest integer lower than m\. The energy spectrum is bounded and again 
not equispaced 

eo < ei < 62 < 63 < . . . < Cm < CAT^ , 

em+i -em = l- ("^+2) ^• 

It can be seen that the energy Cm, considered as a function of m, is an increasing function for low 
values of m and has a maximum at the point 1/A. Beyond that point it becomes a decreasing 
function and it even takes negative values; but, due to the above condition m < m\ for physical 
eigenfunctions, this behaviour of e^ for large values of m is not physically meaningful. Figure III 
shows the form of Cm as a function of m for two particular values of A. 

The following particular cases, that are presented in a decreasing way, illustrate the situation for 
several positive values of A 

• If A > 1 the only bound state of this deformed oscillator is the fundamental state ^o- 

• If 1 > A > 1/2 there are two eigenfunctions; the fundamental level ^0 and another bound 
state ^1 with energies cq < ei. For example, for A = 0.8 we have cq = 1/2 and ei = 1.1. 

• If 1/2 > A > 1/3 there are three eigenfunctions ^O) ^i) ^2) with energies cq < ei < 62- For 
example, for A = 0.4 we have cq = 1/2, ei = 1.3, and 62 = 1.7. 

• If 1/3 > A > 1/4 there are four eigenfunctions ^O) ^i) ^2i ^3) with energies eo < ei < 62 < 
63. For example, for A = 0.3 we have eo = 1/2, ei = 1.35, 62 = 1.90, and 63 = 2.15. 

• Ifl/(n — l)>A>l/n there exist n eigenfunctions (the fundamental one and n — 1 other 
excited bound states) ^0 and ^j, i = 1, 2, . . . , n — 1, with energies cq = 1/2 < ei < 62 < 63 < 
... < e„_i. 

We close this section with a comparison of the energy levels em, (A) of the A-dependent oscillator 
with the corresponding values of the linear harmonic oscillator 

e„(A) = e„(0) + im2 |A| , 



e^(A) = e^(0)-im2A 




(19) 

_ m2 A A ^ n I 
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It is clear, therefore, that when A is negative the energy em(A) is higher than the energy 6^(0) of 
the harmonic oscillator and when A > 0, em(A), m < ttia, is lower than 6^(0). Figure IV shows 
this property by plotting the values of e„j(A) for A = —0.30 and A = 0.30. 



3 Schrodinger factorization formalism 

The factorization of Schrodinger Hamiltonians in terms of first order differential operators [41]- [45] 
has been shown to be very efficient in finding properties of the spectrum of the Hamiltonian and 
even exhaustively solving the spectral problem in the case of shape invariant Hamiltonians [46]- [53]. 

In the usual case one starts with a Hamiltonian Hi like 



H 



dx^ 



+ Viix) 



determined by the potential Vi and looks for two operators a and a~^ of the form 



a=A + »'W, 



|: + '*'(^' 



such that H] admits the factorization 



H^ = a^ a 



d 
dx 



+ W{x) 



^ d_ 
dx 



+ W{x) 



(20) 



This happens if and only if the function W is such that 

i.e., ly is a solution of the preceding Riccati equation. When Hi admits such a factorization we 
can define a new Hamiltonian H2 , of the form 



Ho 



dx"^ 



+ V2{X) 



by means of the new factorization 

H2 = aa^ 



i. + „.,,,] [_^ + „,-(,) 



(21) 



Such a Hamiltonian H2 is said to be the partner of Hi, and leads to a new Riccati equation 

W' + W'^ -V2 = Q. 

In summary, the function 1^ is a solution of two different Riccati equations, determined by the two 
potentials Vi and V2: 

W' -W'^ + Vi = Q, 
W' + W'^ -V2 = Q, 
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Therefore, the two potentials Vi and V2 are related by 

V2 = 2W'^ - Vi . 

The important point is that a is an intertwining operator for both Hamiltonians and this property 
allows to establish relations between the corresponding spectral problems. 

The method can be adapted to deal with the class of Hamiltonians we are interested in: 

cF , d 



Hi 



{l + \x' 



dx"^ 



\x—- + Uiix), 
dx 



which correspond to a problem with a position dependent mass; this is a particular example or 
factorization of Sturm-Liouville operators [59]. In this case we should look for two operators A 
and A^ of the following form 



A 

A+ 



VT+A 



x^ 



■ d_ 
dx 



+ W{x) 



Vl + Ax 



2 — + W{x) , 
dx 



and then, as 



H^ =A+A 



Vl + X: 



dx 



+ Wix) 



vT+x 



x^ 



dx 



+ W{x) 



(22) 



such a factorization is possible if and only if the function W satisfies the Riccati differential equation 



The partner Hamiltonian is now of the form 



Ho 



[l + Xx^ 



dx"^ 



Xx^ + U^ix). 



and is defined by the alternative factorization 



H2 = AA-^ 



d 



Vl + Xx^ — + W{x) 
dx 



d 



Vl + Ax2 — + VF(x) 
dx 



which leads to 



Vl + Xx^W' + W^ -U2 = 0. 



(23) 



In other words, the function W defining the factorization should be a solution of two Riccati 
equations defined by the two potential functions Ui and U2 



Vl + Xx^W' -W^ + Ui = 0, 
Vl + Xx^W' + W^ -U2 = 0, 

and therefore both potentials are related by 

U2 = 2W^ - U^ . 
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4 Schrodinger factorization approach to the A-depen- 
dent nonhnear oscillator 

In a recent paper [40] a factorization for the A-dependent quantum nonlinear oscillator was proposed 
and it was shown to have the shape-invariance property, therefore the spectrum of such system has 
been fully determined. We aim in this section to do a more complete computation of both the 
spectrum and the corresponding eigenfunctions and to study some interesting properties of the 
polynomial replacing Hermite polynomials in this deformed case. 

4.1 Operators A and A'^ and intertwined Hamiltonians Hi and H2 

It has been proved in [40] that the following operators A and A~^ 



A 

A+ 



h 



'1m 

n 

'2m 



d 



,ma. 



ax h y/l + Xx^ 



d 



,ma. 



dx ^ n ' Vl + Ax2 



provide a factorization for the A-dependent quantum nonlinear oscillator. In fact, if we denote by 
W\ the following function 



Vl + Ax2' 



then the product A'^ A that is given by 



A+ A 



2m 



VT+A 



a;2 — + {^-) Wx 
dx n 



Vl + A: 



d ,ma. 



dx 



+ 



h 



Wx 



turns out to be 



A+A 
Ui 



2m 



(1 + Ax' 



2, ^2 d 

2^ + Ax 



(ix2 



dx 



+ Ui 



-{ma)ia + -X)Wi--iha). 



(24) 



In a similar way, the product AA^ given by 



2m 



r r 7; d /W-axTT^ 



r r — 7; d ,fna.^^^ 



turns out to be 



AA^ 
U2 



2m 



(1 + Ax' 



o, (i2 , d 
2^ + Ax 



(ix2 



dx 



+ U2 



{^ma){a--\)Wl + -{na). 
2 m, 2 



(25) 



Consequently, the Hamiltonian H given by (6), is such that 

H = Hi + ]^ (ha) =H2-\ (ha) 
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where Hi and H2 are the Hamiltonian Hi = A^ A and its partner H2 = A A'^ . 
Note that the hmits when A ^ are given by: 

^ h^ (f I 1 

^ h^ (p I 1 

and note also that the commutator [A, A'^] is not a constant. In fact, 

but when A ^ we find 

hm A^o [A, A~^]=ha. 

4.2 The spectrum of the quantum deformed nonUnear oscillator 

The two operators, the Hamiltonian Hi and its partner H2, depend on the parameter a; the 
important point is that this a-dependence determines an important relation between them given 

by 

H2{a) = Hi{ai) + R{ai) , (26) 

where ai denotes ai = /(a) with the functions / and R defined by the expressions 

h 1 h"^ 

f{a) = a A, R{a) = ha + {-)—X. (27) 

m 2 m 

The condition (26) is called 'shape invariance' condition and for such cases there exists a recipe 
to compute the spectrum and to determine the corresponding eigenfunctions [47, 48] (see also [46] 
for a modern approach). 

More specifically, as a first step, the bound state |^o(ct)) is found by solving A{a)\^o{a)) = 0, 
and has a zero energy, i.e. Hi{a)\^o{a)) = 0. Then, using (26) we can see that |^o(c>^i)) is an 
eigenstate of H2{a) with an energy Ei = R{ai), because 

^2(a)|^o(ai)) = (Hiiai) + i?(ai))|^o(ai)) = i?(ai)|^o(ai)) • (28) 

Next, ^^(0)1^0(01)) is an eigenstate of Hi{a) with an energy Ei = R{ai), because if we use the 
intertwining relation Hi[a)A^{a) = A^{a)H2{ct) we see that 

Hi{a)A^a)\^o{ai)) = A\a)H2{a)\^o{ai)) = A\a){Hi{ai) + i?(ai))|^o(«i)) , 

and hence we arrive to 

Hi{a)A\a)\^^{ai)) = R{ai)A\a)\^o{ai)) • 
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This process should be iterated and we will find the sequence of energies for Hi{a) 

k 
Ek = Y.R{aj), Eo = 0, (29) 

i=i 

the corresponding eigenfunctions being 

|*„(ao)) = At(ao)At(ai) • • • A\an^i)\^oian)) , (30) 

where ao = a and Oj = f{aj^i), namely, aj = /■'(ao) = /■'(o^)- 

We can now apply this procedure with the functions / and R given by (27). In this case we 
obtain 

otk = J {a) =a A; , 

m 

R{ak) = hak + ^— \ = n{a k) + ^— \ (31) 

and therefore the energy £"„ that is given by 

becomes 

( fi\ \ nil? { r?'h\\ 

En= \na- --—n{n+ l)\n + - — A = [na 71, 

V Im J Im V Im ' 

that in terms of the adimensional variable A becomes 

En 



n-{-)r? K {ha), n = 0,1,2,... (32) 



^2 
The eigenstates of H will include an additional term {l/2)fia which shifts all levels. 

4.3 Eigenfunctions of the quantum deformed nonlinear oscillator 

As indicated before, the fundamental state is determined by A|^o) = 0. The wave function ^q{x) 
in the coordinate representation is then given by the normalized solution of the linear first order 
equation: 

(l + Ax^) — ^o + /3x^o = 0, /5=^' 

whose solution is proportional to 

1 =A = ± 

°~(1 + Ax2)n)' '^°~2A~2A' 

In the limit of A going to zero, having in mind that \i\iv\-^q{1 + \x^) = 1 and liniA— >o/?/(2 A) = oo, 
we have that 



lim(l + Xx^r^l'' = exp {- lim Ml±A:^') 



13 2 

exp ( --X 



17 



The other eigenfunctions are obtained according to the recipe we derived before as follows: 



*i = A+(ao)^o(ai) 

*2 = ^+(ao)^+(ai)^o(a2) 

^3 = ^+(ao) ^+("1)^^(02) ^0(03) 



*r 



A+{ao) A+{ai) ^+(02) • • • A+(an-i) ^o(an) 



(33) 



In the next property we establish a powerful tool for computing the action of different operators 
A'^{ak) on functions ^o(q^j)- 



Proposition 2 Let g{x) he an arbitrary differentiable function. Then the following equality holds: 



d 



(1 + Ax^f VT+A^-^ (1 + \x^)-Pg{x) 
dx 



■11 



vr+Ax2 — + 



d (2pA) X 



dx Vl + Ax2 



g{x)- 



Proof: This property is proved by direct computation. 
Using the notation Zx = 1 + Ax^, the above property can be rewritten as follows 

^ ^ d 



Zrr, \/ Zx 



dx ^ 



■1 



d ^ {2pX)x 
^ dx ^/% 



in such a way that in the particular case p = b/2X it reduces to 



,6/2A^/^_l^-b/2A 



Zx V^x J ^x ' ~ { ^) 



d bx 



'Zx h 

dx 



The operators A^{ao), A'^{ai), ^+(02), • • ., which are given by 



A+{ao) -- 


= Co[ 


, u 


A+(«i) = 


- C,[ 




A+(«2) = 


- C2[ 


^ d 



,mao x 



' h yTZx 

,ma\ X 



' h JYx 

,ma2. X 



h 



I Zrf 



can be rewritten as follows 



A+(qo) 
^+(ai) 
^+(02) 



C,{-X)zM^^^xl^Zx^^^l^\ 
C^i-l)z^x^/''V^x^z-^^/^\ 
C,{-l)zt"^^xl^Zx^'^/'\ 



where j3q and /3fc denote 13q = (3 and [5k = {mak)/h. 
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Finally, the functions ^O) ^i) ^2, ^3, • • •; ^n, • • •, turn out to be 

*3 = {C,...C,){-Xfz^JI^\zll^.^^^{^fz-J^I^\, 



and therefore the n-th wave- function ^n(2;) is given by 



^„ = (Co...C„)(-l) 






dx^ 



z^.i'-'^z-zM^ 



-P/2\ 



(34) 



that must be considered as a new way of representing the expression (9) obtained in the Sec. II by 
the direct approach to the differential equation. 



5 A-dependent Hermite polynomials 



We have solved the A-dependent nonlinear oscillator by using two different procedures and in both 
cases we have arrived to an expression for the wave- function ^n(x) depending of some A-dependent 
polynomials that must be considered as a deformation of the Hermite polynomials. A natural idea 
is to expect that the classical properties of the Hermite polynomials must remain true, with the 
appropriate modifications, for these new A-dependent Hermite polynomials. In Sec. II we proved 
the orthogonality, now we will study three other fundamental properties: A-dependent Rodrigues 
formula, the existence of a A-dependent generating function and the existence of recursion relations 
among polynomials of different orders. 

The Schrodinger factorization approach to the A-dependent nonlinear oscillator, studied in Sec. 
IV, has provided the expression (34) for the wave function ^„. Therefore, we have 



H„(x,A) = (-l)".fA+i/2 



dx"^ 



or 



H„(y,A) 



l)nU/A+l/2 « 



,n 2- (/3/A+1/2) 



n -(l/K+l/2) 

Zy Zy 



n = 0,l,2,... 



1 + Ay2, 



(35) 



that must be considered as the "Rodrigues formula" for the A-dependent Hermite polynomials 
obtained in Sec. II. In fact, it is clear that 



lim 



and consequently 



A^O 



i\n ^l/A+1/2 _^ 

^ ^ ' dy 



,n ,- (l/A+1/2) 

Zy Zy 



jn 

(-l)"e^^e- 



limA^oWn(y,A) = i:f„(y), n = 0,1,2,... 

In a similar way to the standard Hermite polynomials, every function Tim is a polynomial of degree 
m and it has a well defined parity 



-Hr 



in I m— 2 I m— 4 , 

CmV + Cm-2y + Cm-iV + ■■■ 
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in such a way that the coefficient Cm of x™ in Tim is given by 

Orri 1 

Cm = n - (2-rA) = (2-mA)(2- (m + 1) A) . . . (2 - (2m - 1) A) 

so that 

liniA^o Cm, = 2"" . 

The first polynomials Tin = T~iniy, ^)' f^ = 0, 1, 2, . . . , 6, have the following expressions 

-Ho = 1 
Til = hy 

Ti2 = fc2[2(l-A)y2-l] 

Tis = A;3[2(l-2A)y3_3y] 

Tii = A;4[4(l-2A)(l-3A)/-12(l-2A)|/2+3] 

Ws = fe5[4(l-3A)(l-4A)y5-20(l-3A)/ + 152/] 

Tie = fcg [8(1 - 3A)(1 - 4A)(1 - 5A)/ - 60(1 -3A)(1-4A)/ + 90(1 -3A)y2_ 15] 

where the constants ki, i = 1,2, ... ,6, are given by 



ki 


= (2-A), 


fe2 = (2 - 3A) , 


ks 


= (2-3A)(2-5A), 


fc4 = (2-5A)(2-7A), 


k5 


= (2-5A)(2-7A)(2-9A), 


ke = {2- 7A)(2 - 9A)(2 - llA) . 



We next analyze the existence of a generating function. 

A generating function for the A-dependent Hermite polynomials Tin{y,^ must be a function 
J-{t, y. A) that generates such polynomials as the successive coefficients of its Taylor power series 
of the additional variable t. Since the generating function of the Hermite polynomials is known to 
be given by 

nt, y) = TZ=,^h^ ^«(^) *" ' ^(^' y) = ^^'*'"*'^ ' 

then, a necessary condition to be satisfied by J-{t, y. A) must be the fulfillment of the limit 

limA^o^(i,y,A)=e('*^"*'). 
We propose to consider the function 

T{t,y,A) = (l + A{2ty-t^)Y^^ (36) 



as the appropriate function for representing the A-dependent generating function. It is clear that 
it satisfies the correct limit 

limA^o(l + A(2ty-t2))^/^ = eP*^-*^). 
The power series of J-' in the new variable t is given by 

(l + A (2ty - t') ) '^^ = Er=o(;7! ^ ^"^^' ^^ *" ^^^^ 
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where we have used the notation Tin for the coefficient of the Taylor series. The first polynomials 
Tin = Ti-niy, ^)) n = 0,1,2, ... ,6, obtained in such a way, have the following expressions 

Ho = 1 

Hi = 2giy 

7^2 = 2<72[2(1-A)/-1] 

Hs = 4<73[2(l-2A)y3-3y] 

7^4 = 454[4(l-2A)(l-3A)/-12(l-2A)/ + 3] 

^5 = 8<75[4(l-3A)(l-4A)/-20(l-3A)y3 + i5y] 

■Hq = 856 [8(1 - 3A)(1 - 4A)(1 - 5A)/ - 60(1 -3A)(1-4A)/ + 90(1 -3A)y2- 15] 

where the constants gi, i = 1,2, ... ,6, take the values 

51=52 = 1, 53 = 54 = (1-A), 55 =56 = (1-A)(1-2A). 

Hence, the polynomials Tin, obtained from the generating function J-'{t, y. A), are essentially the 
same ones that the polynomials Tin, obtained from the Rodrigues formula. They coincide in the 
fundamental part (y-dependent polynomial part written between square brackets [...]) and only 
differ in the values of the global multiplicative coefficients that now are given by gi instead of ki 
(although gi and ki are rather similar nevertheless gi ^ ki). 

The existence of this generating function allows us to establish the existence of A-dependent 
recursion relations among deformed Hermite polynomials of different orders. 

Firstly, let us take the derivative of J-^{t, y, A) with respect to t and then multiply both sides by 
1 + A {2ty — t^). In this way we obtain the relation 

that leads to 

Hi = 2y?Yo 
and to the following recursion relation 

?Yn+i = 2y(l-nA)?Y„-n(2-(n-l)A)?Y„_i, n> 1. 

It is clear that if we particularize this formula for A = then we recover the recursion formula of 
the standard Hermite polynomials. 

Finally, if we take the derivative of the function J-{t, y, A) with respect to the variable y and then 
multiply both sides by 1 + A {2ty — t^) we arrive to the following relation 

2t (1 + A {2ty - t')) '"" = (1 + A {2ty - t")) ( £ 1 ^^(y, A) A 

so that, equating the coefficients of each power of t, we obtain 

^0 = . 

H'2 + ^Ayn[ - 2AyHo = 2Ho 

as well as the general expression 

W^+2 + (n + 2)A[2y?Y,l+i-(n+l)H^] = 2(n + 2)?Y„+i , n = 0,l,2,... 
Once more we obtain for A = a well known property of the standard Hermite polynomials. 
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6 Final comments and outlook 

We have have studied the quantum version of the A-dependent one-dimensional non-hnear oscil- 
lator of Mathews and Lakshmanan. The first dificulty to deal with was the question of the order 
ambiguity in the quantization of the Hamiltonian since it is a system with a position-dependent 
mass; this problem was solved by introducing a prescription (previously studied in Ref. [40]) ob- 
tained from the analysis of the properties of the classical system (existence of a Killing vector 
and a A-dependent invariant measure). The quantum Hamiltonian obtained by this procedure has 
proved to belong to the restricted family of exactly solvable Hamiltonians. This is a very inter- 
esting situation: quantum systems with nonpolynomial rational potentials (and standard kinetical 
term) are hard to study and difficult to be solved and, on the other hand, quantum systems with a 
position-dependent mass are also difficult to study and generally not exactly solvable. The system 
we have considered combines however both traits in a very specific way, and this combination turns 
the problem exactly solvable. 

We recall that has been proved in [3] that although the classical two-dimensional nonlinear 
oscillator is not separable in cartesian coordinates it admits however Hamilton-Jacobi separability 
in three other different coordinate systems (polar coordinates and two other A-dependent systems). 
It is known that classical Hamilton-Jacobi separability and quantum Schrodinger separability are 
closely related, so it is to be expected that the n = 2 quantum system will also admit a solvable 
Schrodinger equation. We think that this question must be investigated, as well as its connections 
to other approaches related with the factorization of two-dimensional systems [55], [58], [60]. For 
instance, it was shown in [60] that the Calogero and Calogero-Sutherland models [61, 62] have the 
generalized shape invariance property. 

It is intresting to remark the existence of another family of Hermite-related polynomials connected 
with the harmonic oscillator: the so called "relativistic Hermite polynomials" introduced in [63] 
in the study of quantum relativistic harmonic oscillators (see also [64] and [65] and references 
therein). Although our approach has been entirely non-relativistic and the "A-deformed Hermite" 
polynomials Tin (or Tin) seem to be rather different to these "relativistic Hermite polynomials", 
the existence of some kind of relationship seems quite probable and would be worth studying. 

Finally, also in Ref. [3], a geometrical interpretation of this nonlinear system was proposed 
(a similar interpretation is discussed in [66, 67] with A in the kinetic term and a A-independent 
potential). The main idea is that this A-dependent system represents a nonlinear model for the 
harmonic oscillator in the circle S^ (A < 0) or in the hyperbolic line H^ (A > 0). Remark that 
in the n = 1 case, it is convenient to consider S^ and H^ as one-dimensional spaces obtained by 
endowing each single geodesic of 5^ or H^ with the induced metric (motion on S^ and H^ will 
correspond to the J = radial motions on S'^ or H^). The point is that, if we assume that this 
geometric interpretation is correct, then the A-dependent "deformed Hermite" polynomials Tin (or 
Tin) could be interpreted as the "curved" version of the Hermite polynomials on the circle S^ or in 
the hyperbolic line H^. We think that this is another question to be studied. 
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Figure Captions 

• Figure I. Plot of Vx{x) = (1/2) (a^a;^)/(l + Ax^), a = 1, A < 0, as a function of x, for 
A = — 2 (upper curve), and A = — 1 (lower curve). 

• Figure II. Plot of Vx{x) = (1/2) (a^x^)/(l + Ax^), a = 1, A > 0, as a function of x, for 
A = 1 (upper curve), and A = 2 (lower curve). 

• Figure III. Plot of the energy Cm as a function of m for A = 0.30 (lower curve) and 
A = 0.15 (upper curve). The curves also show the plot of the points (m, Cm) for the values 
m = 0, 1, 2, 3, and m = 0, 1, . . . , 6, respectively. For A = 0.30 there exist four bound states 
(four thick points in the curve) and for A = 0.15 seven bound states (seven thick points in 
the curve) . When A decreases the maximum of the curve moves into the up right and in the 
limit A ^ the curve converges into a straight line parallel to the diagonal (dashed line). 

• Figure IV. Plot of the energy Cm as a function of m for A = 0.30 (lower curve) and 
A = —0.30 (upper curve). The thick points {m, Sm), corresponding to the values m = 0, 1, 2, 3, 
represent the four bound states existing for A = 0.30 and the first four bound states for 
A = —0.30. The straight line (dashed line) placed in the middle corresponds to the linear 
harmonic oscillator. 
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Figure I. Plot ofVxix) = (1/2) (q^x^)/(1 + Ax^), a = 1, A < 0, as a function of x, for A 
(upper curve), and A = — 1 (lower curve). 
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Figure II. Plot of Vx{x) = (1/2) (q^x^)/(1 + Ax^), q = 1, A > 0, as a function of x, for A = 1 
(upper curve), and A = 2 (lower curve). 
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Figure III. Plot of the energy e^ as a function of m for A = 0.30 (lower curve) and A = 0.15 
(upper curve). The curves also show the plot of the points (m, em) for the values m = 0,1, 2, 3, and 
m = 0, 1, . . . , 6, respectively. For A = 0.30 there exist four bound states (four thick points in the 
curve) and for A = 0.15 seven bound states (seven thick points in the curve). When A decreases 
the maximum of the curve moves into the up right and in the limit A ^ the curve converges into 
a straight line parallel to the diagonal (dashed line). 




Figure IV. Plot of the energy e^ as a function of m for A = 0.30 (lower curve) and A = —0.30 
(upper curve). The thick points {m,em), corresponding to the values m = 0, 1,2,3, represent the 
four bound states existing for A = 0.30 and the first four bound states for A = —0.30. The straight 
line (dashed line) placed in the middle corresponds to the linear harmonic oscillator. 
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